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Abstract 

The dynamics of the resonant Auger (RA) process of the core-excited C*0(ls^^7r*, v^. = 0) 
molecule in an intense X-ray laser field is studied theoretically. The theoretical approach includes 
the analogue of the conical intersections of the complex potential energy surfaces of the ground and 
'dressed' resonant states due to intense X-ray pulses, taking into account the decay of the resonance 
and the direct photoionization of the ground state, both populating the same final ionic states 
coherently, as well as the direct photoionization of the resonance state itself. The light-induced 
non-adiabatic effect of the analogue of the conical intersections of the resulting complex potential 
energy surfaces gives rise to strong coupling between the electronic, vibrational and rotational 
degrees of freedom of the diatomic CO molecule. The interplay of the direct photoionization of 
the ground state and of the decay of the resonance increases dramatically with the field intensity. 
The coherent population of a final ionic state via both the direct photoionization and the resonant 
Auger decay channels induces strong interference effects with distinct patterns in the RA electron 
spectra. The individual impact of these physical processes on the total electron yield and on the 
CO^{A ^n) electron spectrum are demonstrated. 
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I. INTRODUCTION 



Recently, the resonant Auger (RA) decay of atoms exposed to strong X-ray laser pulses 
has been studied theoretically ll-4|. It has been demonstrated that under the extreme field 

n n 

conditions provided by X-ray Free Electron Lasers (XFELs) [5|, l6|, such as unprecedented 
high- intensity, ~ 10^^ W/cm^, and very short pulse durations, 10-500 fs, the stimulated 
emission from the resonance back to the ground state starts to compete with the Auger 
decay. The interplay between the resonant excitation and stimulated emission results in 
Rabi fiopping between the ground state and the resonance within its Auger decay lifetime 
and pulse duration. This leads to spectacular modifications of the RA spectra, like the ap- 
pearance of multipeak structure in the electron spectra jl, 2]. It has also been demonstrated 
2|-y], that direct photoionizations of the ground state and of the resonant state increase 
dramatically with the field intensity, resulting in additional leakages of the corresponding 
populations into all possible final ionic states. Moreover, a final ionic state is populated 
coherently by both the direct photoionization from the ground state and by the resonant 
Auger decay inducing strong interference effects with distinct patterns in the RA electron 
spectra of atoms 

What happens when a diatomic molecule is exposed to a strong laser pulse with a carrier 
frequency which fits to the energy difference between two electronic states? It has been 
demonstrated that the two resulting 'dressed' electronic states exhibit a conical intersection 
,n the space of tKe .otat.ca. angle and tKe ,nte..nuc.ea. d.tance fl Q. The e.e.genee 
of this light-induced conical intersection leads to substantial non-adiabatic effects and to a 
strong mixing of rotations and vibrations jsj. Very recently jo], the RA process in molecules 
in intense X-ray laser fields has been studied theoretically. Also here it has been demon- 
strated that an atomic-like picture (fixed internuclear distance and molecular orientation) 
is inapplicable. Not only the vibrational motion plays a role in the decay as in the case of 
weak fields, but the rotational motion becomes very important in strong fields as well. In 
particular, the appearance of the analogue of conical intersections of the now complex (owing 
to a finite life-time of the electronic decaying state) potential energy surfaces (PESs) of the 
ground and 'dressed' resonant states has been predicted in 9|. This phenomenon gives rise 
to strong non-adiabatic coupling and to strong mixing between the electronic, vibrational 
and rotational degrees of freedom of diatomic molecules exposed to strong X-ray laser pulses. 



2 



To demonstrate the non-adiabatic effects, the resonant ionization channel (i.e., excitation 
and decay of the resonance) has been exphcitly included in the numerical calculations of the 
RA process in the HCl molecule studied in Ref. jo], whereas the leakages and interference 
effects mentioned above for atoms were not taken into account explicitly for transparency 
of presentation. 

In the present work we study how the combined action of all the aforementioned compet- 
ing processes evoked by a strong field influence the multidimensional dynamics of the RA 
effect in diatomic molecules exposed to strong X-ray laser pulses. For this purpose we unify 
the previously developed theoretical approaches for RA in atoms j4| and for molecules joj. 
As demonstrated in Ref. j^, the simultaneous effect of the RA decay and of the direct ioniza- 
tion of the involved states results in additional modifications of the Hamiltonian governing 
the nuclear dynamics on the coupled dressed surfaces. One can expect that these modifi- 
cations will even further enhance the non-adiabatic effects caused by the appearance of the 
analogue of a conical intersection of the complex PESs. In our previous study {q], the case 
of excitation of the dissociative HCr(2p~^cr*) resonant state and its subsequent RA decay 
into one of the dissociative final states has been investigated. There, the nuclear dynamics 
accompanying the RA decay result in the fragmentation of of the HCl molecule, and, as a 
consequence, a rather simple RA spectrum appears. The RA decay of the HCl*(2p~^cr*) 
state is however difficult to study experimentally, owing to many overlapping intermediate 
and final electronic states. In the present work we concentrate on the nuclear dynamics 
accompanying the RA decay of a bound resonant electronic state into a bound final ionic 
state of a diatomic molecule, where the corresponding vibrational structures can be resolved 
in the excitation and decay spectra. 

The RA decay of the core-excited C*0(ls~^7r*) molecule is a perfect candidate for our 



purposes. First of all, for weak fields it is one of 



both, experimentally lOMlSj and theoretically 14 



,he most thoroughly studied processes 



2l| . Therefore, all parameters for the 



calculations are available. Second, the weak field excitation spectrum of the core-excited 



resonance consists of three well resolved vibrational states v. 
probabilities of jibout 87%, 12% and 1%, respectively 
= 250 meV 



15 



0, 1 and 2 with excitation 
The energy separation of 



15] between the Vr levels of the C*0 resonance, is about three times larger 



than their natural lifetime width of F 



80 meV 



22l ]. and thus the underlying so called 



lifetime vibrational interference (LVI, 23|) effects are small. The timescale of the nuclear 
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Highly-excited 




FIG. 1: (Color online) Schematic representation of RA effect of the CO molecule exposed to strong 
X-ray pulses (see text for details). 

vibrational motion, t,^ = 2Tx/Uf. ^ 16.5 fs, and that of the RA decay, = 1/Taug ~ 8-2 fs, 
are rather comparable. All of these imply that one can expect measurable fingerprints of the 
nuclear dynamics in the core-excited state when studying strong field effects. Finally, the 
RA electron spectra corresponding to the decay of the resonance to the three lowest X, A, 



and B states of the CO"^ ion do not overlap (the energy separations are about 2 eV |12|. llTj). 
and can thus be studied separately. In the present work we investigate the RA decay into 
the CO^{A ^n) final ionic state which exhibits the richest vibrational structure of all the 
electron spectra in the weak field case [12I, Q • 
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II. THEORY 



The processes relevant to the present study are represented schematically in Fig. [T] The 
strong X-ray laser field with the photon energy u around 287.4 eV couples the ground and 



the resonance states. As demonstrated in Ref. jo], this coupling leads to damped Rabi 
oscillations of the nuclear wave packets between the ground state and the resonance. The 
CO^{A ^n) final ionic state can be coherently populated via the participator Auger decay 
of the resonance and by the direct photoionization from the ground state. The first channel 
is operative at those times when the wave packet is in the resonant state, and the second 
channel - when it is in the ground state. In addition, there are two mechanisms of losses 
relevant to the RA effect in strong fields. The total photoionizations of the ground state 
and of the resonance lead to leakages of the respective wave packets, as indicated in Fig. [1] 
by the vertical helix arrows from the ground state and from the resonance, respectively. 
These processes, except of the Auger decay itself, are operative only when the pulse is on. 
After the pulse expires, the final population of the ground state can be much smaller than 
In the presence of a strong pulse, multiple ionization of a molecule takes place as well 



6|. Different possibilities for that were discussed in Ref. j4|. However, multiple ionization 
processes, in principle, can be experimentally separated from the RA effect as soon as one 
measures energies of ejected electrons and does not intend to measure ions. 

The currently operating XFEL, Linac Coherent Light Source (LCLS), do not produce so 
far a monochromatic radiation and the X-ray pulses consist of many spikes with random 
fluctuations of the frequency, phase and amplitude js, 6|. The impact of these problems on 
the RA effect of atoms has been studied in In addition, the X-ray pulse is reshaped 
during the propagation through a resonant medium {2, sl. In order to reveal the individual 
contributions of all included physical processes, we concentrate here on the physics a single 
molecule undergoes when exposed to a coherent and monochromatic X-ray pulse. For sim- 
plicity we assume a Gaussian-shaped pulse of duration r centered at to {g(t) = e~(*~*")^/^^) 
with linear polarization of the fleld along the z axis: 

S{t) = SQ^t) cos out = Sq g{t) cos out. (1) 

Here Sq is the peak amplitude, and the pulse-shape function g{t) varies slowly on the 
timescale of lixjou. The cycle- averaged intensity of the fleld is given in atomic units via 
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(1 a.u. = 6.43641x10^5 W/cm^) 

m = ^{£o9{t)}\ (2) 
ana 

where a = 1/137.036 is the fine structure constant. For brevity of presentation, the following 
terminology is used throughout the manuscript. A pulse of a peak intensity well below 1 
a.u. (6.43641x10^5 W/cm^) is referred to hereafter as a weak pulse. A weak field does not 
induce noticeable Rabi oscillations (leads to negligibly small stimulated emission) and can 
be treated perturbatively. Pulses with peak intensities comparable to or larger than 1 a.u. 
are referred to hereafter as strong pulses. Strong fields govern photon absorption processes 
in the extremely nonlinear regime, e.g. by Rabi fiopping 

To describe the dynamics of the RA effect theoretically as a function of time, one needs 
to solve the time-dependent Schrodinger equation for the molecule and its interaction with 
the field (atomic units e = mg = h = 1 are used throughout) 

z4/{t) = nmit) = (h^uc + H,i + b8{t)) ^>{t). (3) 



where D represents the dipole transition 



combines those reported in Refs. 



m 



24|. 



operator. The present approach to solve Eq. 



A. Nuclear dynamics Hamiltonian 

We start with the randomly oriented rotating CO molecule in its ground ro-vibronic state 
CO(X ^S"'",fo = 0, Jo = 0), i.e., with an isotropic distribution over the rotational angle. 
We shall demonstrate below that the non-adiabatic effects of coupling the vibrational and 
rotational degrees of freedom of diatomic molecules induced by the laser field will change this 
distribution dramatically. Following Ref. 2J] the total wave function of a rotating molecule 
as a function of time can be represented via the following ansatz including three electronic 
states: the ground state and the resonance with the electronic wave functions $/ and 
respectively, and the final ionic state CO+(A ^n) plus out going electron of energy e with 
the total electronic wave function 

*(t) = \^i{t))<^i + \^R{t))^R + j \^A{e,t))^\de. (4) 

Here |\l//(t)), |\l/R(t)), and |\1/a(£^, t)) are the time-dependent wave packets propagating on the 
PESs of the ground, resonant and final ionic states, respectively. It should be remembered, 
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that these wave packets depend expUcitly on the nuclear coordinates R and 6. For brevity, we 
will explicitly show the dependence on the nuclear coordinates only in the final Hamiltonian 
matrix ([7]) governing the nuclear dynamics. 

In order to obtain the final set of equations for the propagation of the nuclear wave 
packets, we substitute the ansatz (jl]) in the time-dependent Schrodinger equation and 



project the result onto each electronic state. Similarly to Re 



rotating wave approximation 



25| . the local approximation 



26 



J3,|9 



24| . we imply the 



27| . and redefine ('dress') 



the time dependent wave packets of the resonant and final ionic states as follows 



^rW) = |^i?(t))e+^"* and \^A{e,t)) = \^A{e,t))e 



(5) 



Utilizing the rotating wave and the local approximations is well justified for high energy 



resonant processes considered here 



a 



^. Details of the derivation can be found in Refs. 



In order to derive the final Hamiltonian matrix for the nuclear dynamics accompanying 
the studied process, we notice that the final ionic state CO~^{A ^11) can be populated via 
Auger decay of the resonance by the emission of either sa or s6 RA electrons (as indicated 
in Fig. [1]). These channels form the ^11 symmetry of the overall final state of the ion plus 
RA electron which is the symmetry of the resonant state. Contrary to that, the direct 
population of the CO^(A ^11) state from the ground state may proceed via the emission 
of ea, en or e6 photoelectrons (see Fig. [1]). The direct photoionization channels with ea 
and e6 photoelectrons contribute to the same overall final state as the resonant channels, 
which thus superimpose and interfere. Emission of en photoelectrons, however, results in 
the overall final state of ^S"*" symmetry, which cannot be populated resonantly. Moreover, as 
will be shown below, these transitions are operative via different components of the dipole 
transition operator D. That is why it is important to separate the final states explicitly in 
the ansatz (jl]) (the last integral in Eq. (j4]) includes both types of final states). 

Let us collect the individual nuclear wave packets contributing to the total wave function 
(jl]) into a single vector 



|vl/(e,t)) 



l*r^(^,t)) 



(6) 
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where we have spht the contributions to the final ionic state CO~^{A ^11) produced by the 
emission of e{a/6) electrons and those by the emission of en electrons. We define the matrix 



U{R,9,t) = f{R,e) + 

( Vi{R) - ^Vphit) {Dlit) - {W\t)) sine \ 

{D.,{t)-lW{t))sme VR(i?)-|[r„„, + r,(t)]-a; 
4(t)sin^ V VA{R)+e-oj 

V 4(t) cos ^ Va{R) +e-u ) 

(7) 

With these notations, the final set of equations describing propagation of the nuclear wave 
packets ^ now takes on the following compact form 

z|^(£,t)) = H(t)|^(£,t)). (8) 

The matrix H can be viewed as the effective Hamiltonian governing the two-dimensional 
nuclear dynamics in the considered RA decay process in an intense laser field. Below we 
summarize the physical meaning of each term of the Hamiltonian matrix ((Tj) and provide 
their explicit expressions. 

The T(i?, 9) matrix is the common nuclear kinetic energy operator for the vibrational 
motion along the internuclear distance R and the rotational motion described by the angle 
9 between the polarization vector of the laser pulse and the molecular axis. The functions 
Vi{R), Vr(-R), and Va{R) on the diagonal are the potential energies of the corresponding 
electronic states. The transitions between the electronic states are given by the following 
matrix elements of the total Hamiltonian ([3]) 

($^|i^(t)|$,) = (27r|x|2a)^^e-^-* = D,{t)e-'-\ (9a) 

{^f/'\H{t)\^n) = (2a£(a/5)|l/fi2|l7r27r) = V, (9b) 
($f /')|i^(t)|$,) = {e{a/6)\x\ln)^^e-^-' = d,{t)e-'^\ (9c) 

= (e7r|z|l7r)^^^e-'"* = 4(t)e-'"*. (9d) 



In Eqs. ( I9aj) . (I9cll . and ( l9dl) . the rotating wave approximation [25| has already been utilized, 
and, in contrast to the rapidly oscillating factor e"*"^*, the functions D^it), dx{t) and dz{t) 
vary slowly on the timescale of In ju). The matrix element (l9aj) describes the excitation of the 
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resonant state via the x component of the dipole transition operator. The matrix element 
(I9bp represents the Auger decay of the resonance into the final ionic state CO^{A ^11) with 
the emission of e{a/S) electrons, including both direct and exchange Coulomb integrals. The 
matrix elements fl^ and flQdl) correspond to the direct ionization of the ground state into 
the final ionic state CO+{A ^U) with the emission of s{a/6) or ett electrons, respectively, via 
corresponding components of the dipole transition operator. Here and below, all transition 
matrix elements are assumed to vary slowly with the energy across the resonance and with 
the geometry within the Franck-Condon region, and are replaced by their mean values. 

The potential energy of the ground state Vj{R) in the Hamiltonian matrix ([7]) is aug- 
mented by the imaginary time-dependent term —^Tphit) describing the leakage from the 
ground state due to its direct photoionization into all possible final ionic state (referred 
hereafter as '' GS -leakage^ j4|). Its explicit expression was obtained in j4| in the local approx- 
imation Q and reads: 

r,,(t) = 27rj]|<./,(t)r, (10) 
i 

where d?^i^{t) are the dipole transition matrix elements for the direct ionization of the ground 
electronic state into all possible final ionic states numerated by superscript j, and can be 
computed similarly to Eqs. fl9c|) and fl9dp . The total probability for the direct photoionization 
of the ground state (ITOl) is identical to the quantity %h{t) introduced in Ref. 

T,Hit) = %H{t) = </(t)/c^, (11) 

where o"*^* is the total direct photoionization cross section of the ground state at exciting- 
photon energy u, I{t) is the field intensity ([2]), and the quantity I{t)/u stands in Eq. (fTTll 
for the photon flux. 

Similarly, the total leakage from the resonance with 'dressed' potential energy Vr(-R) — to 
is provided by the imaginary part —^\^aug + r*(^)] on the respective diagonal element of 



H. The first term, — 



a 



aug-! 



describes the usual time-independent leakage due to the Auger 



decay RA-leakage^ [4[). In the local approximation (26|, l27|, the total rate for the Auger 



decay is given by Coulomb matrix element similar to Eq. f lQbp as 24 1 



Taug = 271 Y,\V,\\ (12) 

i 

where the summation over index j runs over all possible RA decay channels. The sec- 
ond term, — |r^.(t), is the time-dependent leakage due to the direct photoionization of the 



resonance RD -leakage'' j4|). As a good approximation {2, 3| it can be chosen equal to GS- 
leakage, Tt,{t) = Tphit), since only outer electrons participate in the direct ionization of the 
ground state and the resonance at chosen X-ray photon energy u. 

As was demonstrated in Ref. jj], the usual direct coupling D^{t) between the ground state 
and the resonance through the laser field [l-3| is augmented by an additional time-dependent 
term —^W{t), named Leakage-Induced Complex coupling LIC-coupling') . This term ap- 
pears only if the photoionization from the ground state and Auger decay are simultaneously 



treated as required. Explicitly the LIC-coupling reads {4] 

W{t) = 27rJ2di{t)V^\ (13) 

j 

where the summation over index j runs over all possible final ionic states accessible by both 
channels. This is an indirect coupling which can be interpreted as follows: The photoelectron 
emitted by the ground state is recaptured by the residual ion to produce the resonance state, 
and reversely, the Auger electron can be captured by the residual ion which then becomes 
the neutral atom in its ground state. The ground and 'dressed' resonance states coupling 
(i.e., the matrix elements H12 and H21) is non-hermitian and operative as long as the pulse 
is on. 

We note that the matrix elements of the permanent dipole moment of the CO molecule 
in the ground and resonance states must also be included in the respective diagonals of 
the Hamiltonian matrix ([7]). However, the permanent dipole moment of CO in its ground 
state is small. In addition, these matrix elements are proportional to the rapidly oscillating 
function e"*"^*, and the average effect of the permanent dipole moment on the pulse duration 
timescale is negligibly small. For simplicity, these matrix elements are not included in the 
present calculations. 

The third row of the Hamiltonian matrix ([7]) describes the nuclear dynamics on the 
studied final ionic state produced via the emission of e{a/6) electrons. The ionic state is 
populated by the direct photoionization from the ground state (matrix element dx{t) sin 6) 
and coherently by the Auger decay of the resonance (matrix element V^), all at a given 
kinetic energy e of the emitted electron, and the created wave packet propagates on the 
'dressed' PESs Va{R) + e — u. The same final ionic state can be produced by the direct 
ionization of the ground state via the emission of en electron. This process is governed by 
the matrix element dz{t) cos 6 in the forth row of the Hamiltonian matrix (J?!). The matrix 
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TABLE I: Parameters for photoionization (PI) and resonant Auger (RA) decay of CO utilized in 
the present calculations. The energy of the exciting radiation is chosen to be w = 287.4 eV. 



State / Quantity / Source 


Value 


Ground state la'^2a'^3a'^4a'^l7r^ba'^ X 




Total direct PI cross section 111. 121 


cr*-?* = 0.24 Mb 


Resonance lf7^2(T^3(T^4o-^l7r^5o-^27r^ ^11 




Total RA decay rate [22j 


r„„„ = 80 meV 

dug ■^KJ AiXV^ » 


Oscillator strength [181 


/ = 0.158 a.u. 


Participator state lcr^2cr^3cr^4cr^l7r'^5a^ A 




Partial Ti A dprav rates [1 81 




EG channel 


ri- = 1.07 meV 


eb channel 


T^fg = 3.46 meV 


Partial direct PI cross sections [18] 




ea channel 


a^^'' = 0.0060 Mb 


evr channel 


a^^"" = 0.0105 Mb 


eb channel 


a^f = 0.0185 Mb 



element H42 = indicates that the resonant population via the emission of en electron is 
forbidden. Separating these two different kinds of processes into two different lines of the 
Hamiltonian matrix ensures that the corresponding wave packets accumulate on the PES of 
the CO'^{A ^n) ionic state incoherently (see also discussion above). 

Finally, the nuclear wave packets \'^^^^^\e,t)) and \'^\{e,t)) contain the information on 
the RA electron spectrum for the production of the 00^(^4 ^11) ionic state. The spectrum 
can be computed as the incoherent sum of the norms of the respective wave packets at long 



times 



24| 



aA{e) = lim(vl/^-/^)(e,t)|vl/(f/^)(e,t)) + lim(vl/:^(e,t)|vl/:^(e,t)). (14) 



B. Intersections of complex energy surfaces 

As mentioned in the introduction, all parameters needed for the present calculations 
of the nuclear dynamics governed by the Hamiltonian matrix ([7]), such as the total and 
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partial direct photoionization cross sections and the RA decay rates, as well as the oscillator 
strength for the resonant excitation, can be found in the literature. These parameters are 
collected in Tab. HI The potential energy curves of the ground, the resonant, and chosen final 



ionic states computed in Ref. 18| are depicted in the upper panel of Fig. [2J The vertical 
excitation energy from the initial state to the Vr = vibrational level of the C*0(ls^^7r*) 



resonance is a; = 287.4 eV 



15|. O 



probability which is about 87% 



15 



all levels, this level possesses the maximal excitation 
The energy curve of the resonant state 'dressed' by 



the field of this energy (i.e., Vr{R) — u) is also depicted in the upper panel of Fig. [2l One 
can see, that this curve crosses the curve of the ground state in the Franck-Condon region 
(see also inset in this panel). 

As has been demonstrated in Ref. jsl, crossing of 'dressed' by field states of a diatomic 
molecule results in the intersections of the corresponding two-dimensional complex PESs. 



It is well known 



28 



29| that for the formation of a conical intersection one needs at least 



two nuclear degrees of freedom whose changes affect the electronic wave function. These 
are not available for a free diatomic molecule. However, two dynamical variables R and 6 
enter the Hamiltonian matrix explicitly and the rotation becomes the missing degree of 
freedom to allow for the formation of intersections for diatomics in the presence of a laser 
field. We point out that the rotational degree of freedom is involved in the nuclear dynamics 
only due to the presence of the laser field and only when the pulse is on. Very important, 
the non-adiabatic couplings, i.e., matrix elements of the nuclear momenta along R and 6, 
between the two 'dressed' electronic states are singular at these intersections jsQl. This gives 
rise to dramatic dynamical effects. 

Without the on- and off-diagonal imaginary corrections to the Hamiltonian matrix ([7]), 
the two-dimensional potential energy surfaces of the 'dressed' electronic states exhibit an 
intersection at 6 = (coupling matrix elements H12 and H21 are proportional to sin^). 
The impact of such light-induced conical intersections on the dynamics of a system has been 
recently demonstrated in Ref. jsj. Due to the presence of RA decay width, leakages and LIC- 
coupling in ([7]), the situation becomes more complicated. The two potential energy surfaces 
in R and 6 space obtained by diagonalizing the electronic Hamiltonian H(i?, 9, t) — T(-R, 0) 
in Eq. ([7]) are now complex and generally exhibit two intersecting points at which the real 
as well as the imaginary parts of the two electronic energies become degenerate 30| . This 
analogue of a conical intersection in the continuum has been named doubly intersecting 
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FIG. 2: (Color online) Upper panel: The potential energy curves of the ground X state of 



CO, the resonant state C*0(ls^^7r* ^11), and the final ionic state CO^{A ^H) from Ref 
together with the curve of the resonant state 'dressed' by the photon energy uj = 287.4 eV. The 
Franck-Condon region for the ground state of CO is marked by the vertical dashed lines and the 
position of curve crossing is indicated by vertical doted line (see also inset). Lower panel: When 
the molecule is exposed to the X-ray laser field of peak intensity /q = 3 x 10^^ W/cm^, the ground 
and resonant potentials become doubly intersecting complex surfaces in R and 9 space. Solid and 
dashed lines represent the projections of the real and imaginary seams of the complex intersecting 
surfaces onto the R9 plane. Two solid circles in the inset indicate the complex intersection points, 
where the real surfaces and the imaginary surfaces are simultaneously degenerate (see also Fig. [3]) . 
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FIG. 3: (Color online) Real (upper panel) and imaginary (lower panel) parts of the resulting two- 
dimensional complex potential energy surfaces. The topology of the surfaces is that of doubly 
intersecting complex energy surfaces (DICES). 

complex energy surfaces (DICES). 

In order to illustrate the DICES appearing in the present case we have diagonalized the 
complex electronic Hamiltonian. The 'dressed' potential energy curves shown in the upper 
panel of Fig. [2] and parameters relevant for the RA decay of CO listed in Tab. [Jwere utilized 
for the peak intensity of the laser pulse of Jq = 3 x 10^^ W/cm^. The real and imaginary 
parts of the resulting surfaces are shown in Fig. [3l The projections of the real and imaginary 
seams of the complex intersecting surfaces onto the R9 plane are shown in the lower panel of 
Fig. m The seam of degeneracy of the real parts of the surfaces is shown by the solid line (see 
also inset for an enlarged scale). The imaginary parts of the surfaces have a complementary 
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seam of degeneracy shown by the dashed hne. The complex surfaces intersect at exactly 
two points which are at the edges of the seams (solid circles in the inset). The present 



case fits well to the discussion of DICES considered in Ref. 



30|. Finally, we note that the 



imaginary parts of the surfaces which reflect the decay of the resulting 'dressed' states are 
now strongly dependent on R and 6 o], [sol (see the lower panel of Fig. ^ in contrast to the 
initially constant values entering the Hamiltonian matrix ([7]). 



III. RESULTS AND DISCUSSION 

In order to illustrate the impact of DICES on the dynamics of the RA decay we performed 
three sets of calculations to which we refer as models. The first model is the exact calculation 
including both vibrational and rotational degrees of freedom in the nuclear dynamics (full 
Hamiltonian matrix ((Tj)). We refer to this model as the ^DICES^ model. In the second model, 
only the vibrational motion is allowed in the nuclear dynamics, eliminating thereby DICES 



(hereafter referred to as the 'no DICES^ model 
matrix ([7]), we apply the common approach 



. To exclude rotations from the Hamiltonian 



31 



33| incorporating rotational transitions 



via the selection rule A J = ±1. The approximate Hamiltonian used for the vibrational 
dynamics of diatomics in laser fields can be obtained in this approach by replacing cos 9 and 
sin 6* functions in Eq. ([7]) by their optical transition matrix elements (Fiol cos 6*1100) ~ \ \ 



and (Yi±i| sin6'|yoo) = yl; respectively, where Yi^ are the usual spherical harmonics. In 
the third model (hereafter referred to as the ^Atomic'' model), the whole nuclear dynamics 
is excluded from the calculations. For this purpose, the internuclear distance R is kept fixed 
at the position of the curves crossing (see Fig. [2]) and cos 6* (sin 6*) in Eq. ([7]) are replaced by 
their optical transition matrix elements as in the no DICES model. 

Within all three models, we also study the individual contributions of different leakage 
mechanisms to the RA effect. For this purpose, we performed the full calculations as well 
as calculations in which we take into account only one or two of the leakage mechanisms. 
We shall address the results of these systematic approximations to the three models as: 

• ^Resonant'' - only the resonant channel is accounted for, i.e., all losses due to pho- 
toionization are neglected; 

• ^Direcf - only the direct ionization channel is taken into account, i.e., the resonance 
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state is excluded from the calculations 



• ^Interference^ - both the direct and resonant channels and the interference between 
them are taken into account, but the resonance state itself cannot photoionize; 

• ' Total' - all mechanisms are taken into account including the direct photoionization 
of the resonance. 

The present calculations were performed for Gaussian-shaped laser pulses of durations r = 
4, 8, and 16 fs which are comparable to the RA decay lifetime of the C*0 resonance of 
Td ^ 8.2 fs. To be able to carry out the two-dimensional calculations on the coupled complex 
surfaces we employed an efficient Multi-Configuration Time-Dependent Hartree (MCTDH) 



method 



34| and code [35 1. 



In order to evaluate the LIC-coupling W(t), one has to include the manifold of all final 
ionic states populated coherently via resonant and direct ionization channels (sum over index 
j in Eq. (IT3l)). These individual contributions of different final states to W{t) can be different 
even by sign. In the absence of necessary data on all final ionic states we have estimated 
and utilized in the calculations an upper bound of the true value of W{t). For this purpose 
we assumed that only one direct dipole transition amplitude dx{t) and only one Auger decay 
Coulomb matrix element V enter Eq. f fT3l) . which, however, correspond to the total direct 
photoionization cross section and to the total RA decay probability (Tab. [I]), respectively. 
Our computations shown below illustrate that for CO the impact of the LIC-coupling is 
moderate for the field intensities investigated. 



A. Total electron yield 

The total electron yield as a function of the X-ray peak intensity calculated in the differ- 
ent models and approximations discussed above is depicted in Figs. HI [5l and [61 In Fig. HJ 
the impact of DICES and nuclear dynamics in the presence of different leakage mechanisms 
is illustrated for a pulse duration of 8 fs. Fig. [5] illustrates the impact of different leakage 
mechanisms within the DICES model for the same pulse duration. The total electron yield 
computed exactly (i.e., within the DICES model without approximations {Total)) is com- 
pared in Fig. |6] for different pulse durations. The total electron yield was computed as in 
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^ using 

lim [ dsA'^Ae,,t)\<ilj{e,,t)) = 1 - lim(vl/,(t)|xl/,(t)), (15) 

where index j in the left-hand side of Eq. (1T5|) must run over all possible final ionic states. 

The atomic-like case, where the nuclear dynamics in R and 6 space is excluded completely, 
is well understood jll, 0]. In the Resonant approximation within the Atomic model, only 
the Auger decay of the resonance and the Rabi oscillations between the resonance and the 
ground state are competing with each other. If the duration of the pulse is comparable or 
shorter than the Auger decay lifetime of the resonance and its intensity is low, the system 
has a finite probability of staying neutral after the laser pulse is over (dashed-dotted curve in 
the uppermost panel of Fig. H] intensities below 10^^ W/cm^). At larger intensities the field 
manages to transfer the whole population from the ground state into the resonance, and 
the atomic-like system can be completely ionized with the probability of 1. This situation 
corresponds to the half-completed Rabi cycles and to maxima in the total electron yield 
arriving at unity jlj] (dashed-dotted curve in the uppermost panel of Fig. H]). At certain 
intensities the atomic system manages to complete several Rabi cycles during the pulse (to 
transfer the population back to the ground state) and the ionization probability drops again. 
These intensities correspond to the minima in the total yield. 

In the Resonance approximation for a molecule, the nuclear motion is additionally in- 
volved in the dynamics of the RA decay. As a result, the oscillations in the total electron 
yield possess much smaller amplitudes. Even in the no DICES model, where only R is the 
dynamical variable and 6 is excluded from the dynamics, the difference from the Atomic 
model is substantial {q] (cf dash-dotted and dashed curves in the uppermost panel of Fig. H]). 
The maxima in the total electron yield are not necessarily touch unity (dashed curve). This 
is because a part of the energy of the absorbed photons is now transferred from the electronic 
degrees of freedom to the nuclear vibrational motion. The deviation from the atomic case 
is even more dramatic if we compare with the results of the full dynamics involving R and 
6 motions {DICES model). In the Resonant approximation within this model, the com- 
plexity of the energy surfaces stems from the RA-leakage mechanism (Auger decay width) 
alone, and the appearing DICES are somewhat different {q] from the full case illustrated 
in the lower panel of Fig. El The seams of degeneracy of the real and imaginary parts of 
the surfaces belong to the straight line in the R6 plane, but not to the contour (see also 
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FIG. 4: (Color online) Total electron yield after exposure of CO to a coherent Gaussian-shaped 
pulse of duration 8 fs computed within different models and approximations discussed in the 
text. These include the atomic-like model (Atomic; dash-dotted curves), rotations free model (no 
DICES] dashed curves), and the exact model accounting for the vibrational and rotational degrees 
of freedom (DICES; solid curves). Upper panel: contribution of only the resonant channel (Reso- 
nant). Middle panel: contribution of both the direct and resonant channels and the interference 
between them (Interference). Lower panel: the contributions of all mechanisms including the direct 
photoionization of the resonance (Total). 
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FIG. 5: (Color online) Total electron yield after exposure of CO to a coherent Gaussian-shaped 
pulse of duration 8 fs computed within the exact model, accounting for the vibrational and rota- 
tional degrees of freedom {DICES model). Shown are results of different approximations discussed 
in the text. These include the contribution of only the resonant channel (Resonant), of only the di- 
rect ionization channel (Direct), and of both the direct and resonant channels and the interference 
between them (Interference), as well as the contributions of all mechanisms including the direct 
photoionization of the resonance (Total). The vertical dotted lines indicate the peak intensities 
chosen for the calculation of the Auger spectra in Figs. [7]and[8j 



the original Ref. 30|). The strong non-adiabatic effects produced by the DICES manage 
to transfer energy from the vibrational to the rotational degree of freedom of a molecule. 
Consequently, the oscillations in the total yield are strongly suppressed from below and from 
above, and the yield is always well below 1 (solid curve in the uppermost panel of Fig. H]). 
We emphasize that the strong impact of rotations illustrated in Fig. H] is not due to the 
rotational motion alone, but rather owing to the light-induced conical intersection (more 
precisely, DICES in the present case) which couples the vibrational and rotational motions 
7H9|. In both Atomic and no DICES models, the coupling matrix elements H12 and H21 in 
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FIG. 6: (Color online) Total electron yield after exposure of CO to coherent Gaussian-shaped pulses 
of different durations (4, 8, and 16 fs) computed within the DICES model which accounts for the 
vibrational and rotational dynamics, employing the full formalism {Total: no approximation). 



Eq. ([7]) carry the factor y | after averaging over the rotational motion (i.e. by the optical 
transition matrix elements of sin^). Therefore, periods of the oscillations in the total yield 
computed in these two models cover larger intervals of the field strength than in the DICES 
model. 

The total electron yield computed in the Direct approximation within the DICES model 
(blue curve in Fig. [5]) illustrates the individual contribution of the GS-leakage mechanism to 
the ionization process of the molecule. The total cross section for the ionization of the outer 
shells of CO 

i'^ph — 0-24 Mb, Tab. [T]) is much smaller than the probability for the ionization 
via the resonant channel. Thus, the total electron yield saturates to 1 due to the GS-leakage 
at much larger intensities (at about 10^^ W/cm^) than in the Resonant approximation. After 
its saturation, the direct ionization channel becomes comparable with the resonant one 4]- 
In the Interference approximation within the DICES model (middle panel of Fig. HI 
and also green line in Fig. [5]), four competitive mechanisms are present in the ionization of 
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the molecule. These are the RA- leakage and GS- leakage discussed above, the interference 
between these two channels, and, of course, the wave packet nuclear dynamics on the strongly 
coupled PESs. The imaginary parts of the resulting complex intersecting energy surfaces 
are responsible for the leakages of the propagating wave packets. As is evident from the 
lower panel of Fig. El the resulting complex energy surfaces are strongly dependent on 
the dynamical variables and on all leakages. Consequently, a strict separation between 
the ionization mechanisms is impossible in the adiabatic picture. A simplified picture can, 
however, be suggested in the diabatic representation (before diagonalization of the electronic 
part of the Hamiltonian matrix ([7])). The nuclear wave packet propagates on the coupled 
complex surfaces of the ground and 'dressed' resonant states. The parts of the wave packet, 
propagating on the resonant surface, decay into the final ionic state via the RA-leakage 
mechanism. The parts of the wave packet, which propagate on the ground states surface, 
decay via the GS-leakage mechanism, i.e., direct photoionization. The competition between 
these two leakage mechanisms results in distinct modifications of the yield in all three models 
(middle panel of Fig. Hj). The oscillations become less pronounced and saturate along the 
trend imposed by the Direct approximation (cf also green and red curves in Fig. [5]). In the 
Atomic model, not all of the oscillations arrive now at unity yield (dash-dotted curve in the 
middle panel of Fig. |4]). This weak effect is a consequence of the complex coupling term 
(LIC-coupling) j4|. The dramatic differences between the computed total yields shown in 
the middle panel of Fig. H] illustrate the impact of the strongly coupled ro- vibrational motion 
via the DICES underlying the processes included in the Interference approximation. The 
differences between the dash-dotted and dashed curves are due to the vibrational motion, 
and those between the dashed and solid curves - due to the rotational motion and DICES. 

The photoionization of the resonance included in the RD-leakage mechanism enhances 
the ionization of the molecule ( Total calculations in the lowermost panel of Fig. H) and also 
black curve in Fig. [5]). It results in further losses of the wave packet propagating on the 
resonant surface, reducing its norm by the direct ionization of the resonance into preferably 
highly-excited final ionic states. The oscillations in the total electron yield now become even 
less pronounced and saturate much earlier as a function of the laser intensity (cf green and 
black curves in Fig. [5]). We note that depending on the pulse duration and the parameters 
of the system, there is a certain interval of intensities where both the GS-leakage and RD- 
leakage mechanisms are comparatively weak Q]. For the CO molecule exposed to a Gaussian 
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pulse of 8 fs, these are the intensities up to about 10^^ W/cm^ (see Fig. [5]). The total electron 
yields computed exactly (i.e., in the DICES model and without approximations (Total)) for 
the three pulse durations of 4, 8, and 16 fs are compared in Fig. O For the shorter pulse of 
4 fs, all leakages are weak up to about an intensity of 2 x 10^® W/cm^, and for the longer 
pulse of 16 fs up to about an intensity of 5 x 10^^ W/cm^ which are twice larger and twice 
lower than the corresponding value of the 8 fs pulse, respectively. The dramatic effects of 
the DICES on the ro-vibrational dynamics persists also in the full calculations (cf, all curves 
in the lowermost panel of Fig. |4]). The collective impact of all leakage mechanisms and of 
DICES results in a strong suppression of the oscillations of the total electron yield and in 
the saturation of the yield as a function of intensity, as illustrated in Fig. \5\ by the black 
curve. 



B. Resonant Auger spectrum 

The RA electron spectra of the CO~^{A ^11) final ionic state at different peak intensities 
of the Gaussian-shaped pulse of 8 fs are depicted in Figs. [7] and [H and those for the shorter 
(4 fs) and longer (16 fs) pulses in Figs. [9] and [10], respectively. Shown in these figures are the 
results of calculations performed within the two models [no DICES and DICES) without 
further approximations (Total), i.e., with all leakage mechanisms and interference effects 
included. The experimental spectrum of Ref. measured in the weak field utilized by 
synchrotron radiation is shown for reference at the bottom of Fig. [7] by open circles. One 
can see from the lower panel of Fig. [7] that in a weak field (see the two spectra computed 
within the no DICES and DICES models and peak intensity of 2.3 x 10^^ W/cm^): (i) the 
structure of the spectrum is independent of whether we have included DICES or not; (ii) 
the contribution of all leakage mechanisms due to direct photoionization is negligible; and 



(iii) the present theory well reproduces the experiment 17|. 

Let us first discuss the effect of DICES on the RA electron spectra in the presence of 
all leakages, as well as the trends in the RA spectrum induced by the field intensity and 
pulse duration. These are illustrated in Figs. [71-[T0l At field intensities beyond the first 
maximum in the total electron yield (Fig. [5]), the stimulated emission from the resonance 
starts to compete with the RA decay, resulting in significant modifications of the electron 
line profiles. In atoms the respective modifications consist in bifurcations of the RA electron 
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FIG. 7: (Color online) The RA electron spectra of CO for the Gaussian-shaped pulse of 8 fs 
duration computed employing the full formalism ( Total) within the two models: no DICES (broken 
curves) and DICES (sohd curves). The peak intensities 2.3 x 10^^ W/cm^, 1 x 10^^ W/cm^, 
2.6 X 10^^ W/cm^, and 1 x 10^^ W/cm^ chosen for the calculations are marked by the vertical 
dotted lines in Fig. [5j The weak field experimental spectrum taken from Ref. {l?! is shown for 
reference by open circles at the bottom. The electron energies corresponding to the energy positions 
of the v' vibrational levels of the CO~^{A ^11) state are marked by vertical dotted lines. 
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FIG. 8: (Color online) The RA electron spectra of CO for the Gaussian-shaped pulse of 8 fs duration 
computed employing the full formalism {Total) within the two models: no DICES (broken curves) 
and DICES (solid curves). The peak intensities are 3 x 10^^ W/cm^, and 1 x 10^^ W/cm^ (see 
also notations in Fig. [7]) . 

peaks [l, 4|. However, the situation in molecules is more complicated than in atoms due to 
the presence of many vibrational and rotational levels of each electronic state involved in the 
process. At an intensity of 1 x 10^^ W/cm^ (lowermost spectra in the upper panel of Fig. [7]), 
the local minima in the line shapes appear at energies where maxima of the intensity were 
observed in the weak field experiments, i.e., around the electron energies corresponding to 
the energy positions of the v' vibrational levels of the final ionic state marked by the vertical 
dotted lines in the figure. At a peak intensity of 2.6 x 10^^ W/cm^ (middle spectra in the 
upper panel of Fig. [7]), the line shapes bifurcate again and their maxima again appear at 
their weak field energies. However, relatively weak local maxima are now present in between 
the large maxima. Up to the peak intensity of 1 x lO'^^ W/cm^ shown in Fig. [TJ the effect of 
DICES is strongly visible, but there are no truly dramatic effects to be seen when comparing 
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FIG. 9: (Color online) The RA electron spectra of CO for the Gaussian-shaped pulse of 4 fs duration 
computed employing the full formalism (Total) within the two models: no DICES (broken curves) 
and DICES (sohd curves). The peak intensities are 2 x 10^^ W/cm^ and 1 x 10^'' W/cm^. 

the no DICES and DICES models. The effect of DICES mainly consists in the broadening 
of the vibrational peaks and in the lowering of the intensities of the maxima in the spectrum 
(cf solid and broken curves in Fig. [7] for all intensities). 

At higher peak intensities (Figs. I814T0] for the 8 fs, 4 fs and 16 fs pulses, respectively), a 
manifold of bifurcations of the line profiles has occurred. In addition, all leakage mechanisms 
evoked by the strong field contribute significantly to the computed spectra. As a result, the 
spectra are very different from the experimental one measured in weak fields . The spec- 
tra cover much broader ranges of the electron energies. Very noteworthy is the appearance 
of the two-hump structure at the low and high energy sides at the largest field intensities 
considered here. It is produced by the strong direct ionization from the ground state, the 
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FIG. 10: (Color online) The RA electron spectra of CO for the Gaussian-shaped pulse of 16 fs 
duration computed employing the full formalism ( Total) within the two models: no DICES (broken 
curves) and DICES (solid curves). The peak intensities are 4 x 10^^ W/cm^ and 3 x 10^^ W/cm^. 

influence of which grows rapidly with the fleld intensity j4| (see also the discussion around 
Fig. [11] at the end of this subsection). Obviously, the spectra exhibit many new features. 
Dramatic difference between the spectra computed within the no DICES and DICES mod- 
els can be seen in Figs. [5HTU1 The spectra computed with and without accounting for the 
rotational motion (or with and without DICES) differ qualitatively in their structures, pos- 
sessing different numbers of resonant features at different energy positions, dissimilar widths 
of peaks and even different relative intensity distributions among them. In some cases (see, 
e.g., the lower spectra in Figs. 18) and [TOj) . the positions of local minima and maxima in the 
spectrum are even alternated. 

Let us now discuss the individual contributions of the various leakage mechanisms to the 
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exact RA electron spectra. Fig. [TT] depicts the spectra computed within the DICES model 
for the 8 fs pulse at the peak intensity of 3 x 10^^ W/cm^. The exact spectrum (Total), 
of course, coincides with that shown in Fig. |H1 The spectrum computed in the Resonance 
approximation (i.e., in the presence of only RA-decay) is depicted by the dash-dotted curve. 
As discussed above, at this field intensity several bifurcations of each vibrational peak in 
the spectrum take place, and the resonant structure of the spectrum is wealthier than that 
observed in weak fields. The individual contribution of the direct photoionization from the 
ground state (GS-leakage) is shown in the figure by the dotted curve. This contribution 
is enormous and the curve is shown on a strongly suppressed scale. Since there is no 
coupling between the ground and the final ionic states (see Eq. ([7])), the line profile computed 
in the Direct approximation remains unchanged with increasing laser intensity and hence 
displays no bifurcations. As demonstrated in Ref. j4|, the intensity of the Direct spectrum 
increases rapidly with the field intensity (see also Fig. |5]). If the direct photoionization 
and resonant spectra were independent of each other, the direct ionization would be the 
dominant mechanism for the population of the CO'^{A ^U) final ionic state at this peak 
intensity (cf dash-dotted and dotted curves in Fig. [TT]and note the 1/10 factor). 

Obviously, since the GS-leakage intensifies with increasing field intensity, the role of the 
underlying interference effects in the spectrum also grows Q]. At the peak intensity of 
3 X 10^^ W/cm^, the interference between the resonant and direct channels [Interference 
approximation; dashed curve in Fig. ITT]) results in a surprisingly strong suppression of the 
direct ionization by the competing resonant pathway. The resulting distinct peak patterns 
in the RA spectrum possess strong asymmetry as seen by comparing the dash-dotted and 
dashed curves in the figure. With the present parameters utilized in the calculations (both 
the resonant and the direct electronic amplitudes have the same signs), the interference is 
more destructive on the low electron energy side than on the high energy side. At larger 
intensities, the interference determines not only the shape of the computed RA spectra, but 
their integral intensities as well (see, e.g., the upper spectra in Figs. [Eland [9]). However, 
as the total electron yield saturates as a consequence of the GS-leakage mechanism, the 
impact of the interference effects saturates too We note that strong interference effects 
persist also in the full (Total) calculations (solid curve in Fig. [TT]) . where the losses of the 
resonance due to photoionization is additionally taken into account. Depopulation of the 
resonance due to its ionization into other usually high-lying final ionic states results in an 
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FIG. 11: (Color online) The RA electron spectra of CO for an 8 fs Gaussian-shaped pulse of 
3x 10^^ W/cm^ peak intensity computed within the DICES model (i.e., including rotations). Shown 
are the results of different approximations discussed in the text. These include the contribution 
of only the resonant channel (Resonant; dash-dotted curve), of only the direct ionization channel 
(Direct; dotted curve), and of both the direct and resonant channels and the interference between 
them (Interference; dashed curve), as well as the contributions of all mechanisms including the 
direct photoionization of the resonance (Total; solid curve). Note that the dotted curve is shown 
on a suppressed scale (factor 1/10) compared to the other curves. 

overall suppression of the spectrum without substantially changing its shape (cf dashed and 
solid curves in the figure). At even higher intensity, however, the peaks seen in the middle 
of the spectrum will decrease and finally may practically become invisible as a consequence 
of the strong photoionization of the resonance state. 

C. Impact of DICES on vibrations and rotations 

Before the pulse arrived, the CO molecules have been in their ground ro-vibronic state 
CO(X ^S"*", vq = 0, Jo = 0) which implies that they were isotropically oriented. Fig. [T^and 
[13] illustrate the change in the vibrational and rotational populations of the neutral molecules 
remaining in their ground electronic state after the 8 fs pulse has expired. In the no DICES 
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model, where the rotational motion is excluded, the pulse couples only the electronic and 
vibrational degrees of freedom and the computed distribution of the neutral CO molecules 
over the rotational angle 9 remains isotropic represented by only the Jq = state. However, 
a part of the energy of absorbed photons is transferred from the electronic excitation to the 
vibrational ones. As one can see from Fig. [T2] {no DICES model; broken curves), at the 
relatively low intensities below Jq = 10^^ W/cm^ the remaining neutral CO molecules stay 
mainly in the ground vibrational state = 0. As the field increases, the excitation of the 
f = 1 vibrational state is possible. At the peak intensity of around 2.5 x 10^^ W/cm^, which 
correspond to the first maximum in the total electron yield computed in the same model 
(lower panel of Fig. H]), the populations of the Vq = and 1 vibrational levels become even 
comparable as can be seen by inspecting the broken curves in Fig. [T21 At higher intensities, 
these populations oscillate: the maxima in the population of the vq = 1 state correspond to 
the maxima in the total electron yield depicted in the lower panel of Fig. |4] by the dashed 
curve. We notice that at the considered peak intensities, the population of the vq = 2 
vibrational level is always below 1% and, therefore, is not shown in the figure. 

Fig. Us] illustrates the changes in the rotational distributions of the neutral molecules 
remaining in the vq = and 1 vibrational levels of the ground electronic state after the pulse 
has expired. The impact of the light-induced non-adiabatic effects (DICES), which couples 
the vibrational and rotational degrees of freedom of the molecule, gives rise to the exchange 
of a very large number of angular momentum quanta during the ensuing dynamics. As a 
result, a wide range of rotational quantum numbers Jq up to about Jq = 30 is populated 
in the f o = and 1 vibrational levels of the ground electronic state. The dipole transition 
selection rule of A J = ±1 for each single absorbed/emitted photon ensures the population 
of only even Jq ro-vibronic levels of the neutral molecules remaining in the ground electronic 
state after the pulse has expired. Fig. [13] shows that higher and higher rotational states 
Jo are excited with increasing field strength, and the populations of the rotational states 
possess oscillations with the peak intensity similar to the total electron yield (solid curve in 
the lower panel of Fig. H]). At the relatively low peak intensities below Jq = 10^^ W/cm^, 
only a few rotational states of small Jq are excited. At higher peak intensities of around 
Jq = 10^^ W/cm^ the populations of low and high Jq quantum numbers become comparable. 
The excitation of higher rotational states Jq results in a strongly anisotropic distribution of 
the neutral molecules over the rotational angle 6 after the pulse is over. We remind that in 
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FIG. 12: (Color online) Impact of the DICES on the vibrational population of the neutral CO 
molecules remaining in its ground electronic state after the Gaussian-shaped pulse of 8 fs duration 
has expired. Shown are the populations of the t;o = and 1 vibrational levels as functions of the 
peak intensity, computed within the no DICES (broken curves) and DICES (solid curves) models 
(Total: no approximations). At the considered peak intensities, the population of the vq = 2 
vibrational level (not shown in the figure) is always well below 1%. The sum of the populations of 
all vibrational levels is normalized to 100% at each peak intensity. 



the no DICES model, the distribution over the rotational angle 6 for the ground electronic 
state remains isotropic as only Jq = appears. The coupling between the vibrational and 
rotational degrees of freedom via DICES changes also the vibrational populations of the 
residual neutral molecules, as demonstrated by the solid curves in Fig. [121 The oscillations 
in the vibrational populations become less pronounced (cf solid and broken curves for each 
vibrational state Vq in the figure). The maxima in the population of the vq = 1 state 
correspond now to the maxima in the total electron yield computed in the same model 
(solid curve in the lower panel of Fig. H]). 

The impact of the DICES on the nuclear dynamics is also documented in the distributions 
of ions jsl, e.g., in the final CO^{A ^H) electronic state of interest (not shown in the paper). 
This effect is, however, difficult to measure owing to the following reasons. First of all, these 
rotational distributions of the final ionic states are different for each energy of the emitted 
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FIG. 13: (Color online) Impact of the DICES on the rotational population of the neutral CO 
molecules remaining in its ground electronic state after the Gaussian-shaped pulse of 8 fs duration 
has expired. Shown are the populations of the Jq rotational levels for the vq = (upper panel) 
and Vq = 1 (lower panel) vibrational levels as a function of the peak intensity computed exactly 
(i.e., DICES model and Total). For each vibrational level vq, the sum of the populations of all 
rotational levels is normalized to 100% at each peak intensity. 

RA electron, and coincidence techniques must be applied. In addition, the ions resulting 
from the RA decay can be further ionized by the strong pulse, and, as a consequence, taken 
away from the coincidence signal. 
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IV. CONCLUSIONS 



The resonant Auger decay of the core-excited C*0(ls~^7r*, = 0) molecule into the 
final ionic state 00^(^4 ^11) in an intense X-ray laser field is studied theoretically. The 
mechanisms responsible for the dynamics of the process are incorporated in the calculations. 
In particular, we investigated the impact of the light-induced non-adiabatic coupling between 
the vibrational and rotational degrees of freedom of a diatomic molecule caused by the 
analogue of the conical intersections of the complex energy surfaces of the ground and 
'dressed' resonant electronic states (briefly, by the light-induced DICES). The non-adiabatic 
nuclear wave packet dynamics accompanying the process competes with different ionization 
mechanisms of the molecule evoked by the laser field: direct ionizations of the ground state 
(GS-leakage) and of the resonance (RD- leakage), and the Auger decay of the resonance 
(RA- leakage). 

The dynamics of the RA process of a diatomic molecule exposed to strong X-ray pulses 
is completely different from that predicted for atomic systems, due to the presence of the 
nuclear degrees of freedom. In the adiabatic representation, the two-dimensional nuclear 
wave packet propagates on the complex (owing to the various leakages) potential energy 
surfaces of the ground and 'dressed' resonant states. Being coupled by the non-hermitian 
interaction induced by the laser field and by the RA decay, these complex surfaces become 
doubly intersecting in R and 9 space (DICES; doubly intersecting complex energy surfaces). 
The non-adiabatic couplings at these intersection points are singular, resulting in strong cou- 
pling between the electronic, vibrational and rotational degrees of freedom. The imaginary 
parts of the resulting intersecting complex surfaces, responsible for leakages of the nuclear 
wave packets from the surfaces, are strongly dependent on both R and 9. As a result, one 
cannot anymore make a clear separation between the involved leakage mechanisms. Finally, 
we remind that the finite pulse duration enforces an explicitly time dependent picture of the 
process. 

The individual impacts of the different underlying mechanisms on the dynamics of the 
process are illustrated on observable quantities like the total electron yield, the RA electron 
spectra, and the vibrational and rotational distributions of the neutral molecules remaining 
in the ground electronic state after the pulse is over. Interference effects and the effects of the 
competition between the different leakage mechanisms are identified and discussed. These 
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leakage and interference effects increase with the field strength. The non-adiabatic effects 
result in dramatic changes of the computed total electron yield and electron spectra. These 
changes are not only of quantitative, but also of unexpectedly strong qualitative nature, 
and can be verified experimentally. Finally, although the initial rotational distribution has 
been chosen to be isotropic, the energy exchange between the vibrational and the rotational 
degrees of freedom via the light-induced DICES results in strongly non- homogeneous angular 
distributions of the neutral molecules surviving the pulse and of the ions produced. In the 
absence of the non-adiabatic effects these distributions would remain isotropic as initially 
prepared. 

The controllability of the light-induced DICES, i.e., of the location of the intersection by 
laser frequency and of the strength of the interstate coupling by the field intensity, makes 
the investigation in intense fields of the resonance Auger process and of decay processes in 
general a challenging and promising new area of research. 
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